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We study how topological phases evolve in the Kane-Mele-Kondo lattice at finite temperatures 
and obtain the topological Doniach phase diagram. In particular, we find an intriguing crossover 
behavior induced by the interplay between the topological structure and electron correlations; the 
topological properties are restored by temperature effects. This restoration can be observed in the 
behavior of the bulk as well as the edge. In the bulk, we observe an increase of the spin-Hall 
conductivity at finite temperatures, while it is zero in the low temperature region. At the edge, we 
observe gapless edge modes emerging with increasing temperature. 


I. INTRODUCTION 

After the discovery of topological insulators, nontriv¬ 
ial properties of topological phases have been extensively 
studied. In these systems, various remarkable observa¬ 
tions have been made because of the presence of gap¬ 
less edge modes arising from topologically nontrivial bulk 
propertieSfii^ These modes induce the quantized Hall 
conductivity in integer quantum Hall systems and topo¬ 
logical magnetoelectric effects in three-dimensional topo¬ 
logical insulators. Furthermore, the edge mode in super¬ 
conductors is the origin of exotic particles. Gapless edge 
modes in the one-dimensional topological superconductor 
are composed of Majorana fermions whose anti-particles 
are themselves, and their experimental realization has 
been addressed recentlyi^^— 

While topological insulators in free-fermion systems 
have been extensively studied, the impact of correlations 
on these systems becomes recently one of the important 
issues in this field^^— because several d- and /-electron 
systems have been proposed as correlated topological 
insulators.ii2r22 In addition, the combination of topolog¬ 
ical structures and electron correlations is expected to 
trigger off exotic phenomena. An example is the topo¬ 
logical Mott insulator—, where nontrivial properties of 
the single-particle Green’s function lead to gapless edge 
modes only in the spin excitation spectrum— 

The topological Kondo insulator— is another exam¬ 
ple of strongly correlated topological insulators. This 
phase attracts much interest recently because there are 
various proposals for topological insulators in heavy- 
fermion systems. It has been proposed that SmBg 
is a realization of such a strong topological insulator, 
where in-gap states have been observed via ARPES 
measurements— YbBi 2 is proposed as a candidate 
material for topological crystalline Kondo insulators— 
The topological phases in heavy-fermion systems have 
been extensively studied at zero temperature) ^°i^^~ — be¬ 
cause the topological structure is well-defined in this 
limit. 

An open question is, however, how finite tempera¬ 
tures affect topological Kondo insulators. This question 


is of particular interest because experiments for topo¬ 
logical insulators are carried out at finite temperatures. 
For example, the Doniach phase diagram, which is well- 
known as a qualitative phase diagram for ordinary heavy- 
fermion systems, has not been established for topologi¬ 
cally nontrivial systems. Furthermore, at finite tempera¬ 
tures, we can expect peculiar behaviors due to the inter¬ 
play between electron correlations and topologically non¬ 
trivial properties because characteristic temperature de¬ 
pendence has been observed even in the absence of topo¬ 
logical properties. Although several groups have studied 
topological Kondo insulators at finite temperatures)^^— 
temperature effects on topological phases have not been 
systematically explored. 


In this paper, we elucidate how topologically nontriv¬ 
ial/trivial phases evolve at finite temperatures. Specif¬ 
ically, we study the Kane-Mele Kondo lattice at fi¬ 
nite temperatures with the dynamical mean field theory 
(DMFT). The obtained Doniach phase diagram shows 
three phases at zero temperatur o^^'^^ : an antiferromag¬ 
netic topological insulator, an antiferromagnetic trivial 
insulator, and a trivial Kondo insulator. Furthermore, we 
find a restoration of topologically nontrivial properties at 
finite temperatures because of the interplay between the 
topologically nontrivial structure and the electron cor¬ 
relations. In order to study these phases, we analyze 
bulk and edge properties. In the bulk, the spin-Hall con¬ 
ductivity which is almost zero around zero temperature 
increases with increasing temperature. At the edge, the 
gapless edge modes emerge with increasing temperature. 
We elucidate that the interplay of the topological nature 
of the system and the Kondo effect is essential for the 
restoration of topological properties. 


The rest of this paper is organized as follows. In the 
next section (Sec. |IT| we explain our model and method. 
The details of the results are discussed in Sec. uni which 
is followed by a short summary in Sec. IIVI 



2 


II. MODEL AND METHOD 

In order to study topological properties in heavy- 
fermion systems, we analyze the Kane-Mele Kondo 
latticei^i^ The Hamiltonian reads 

Hkmk = Hkm + J ^ ■ St, (la) 

i 

Hkm = -^ X] 

{i,3)cr 

+it,o X sgn(i, j)sgn(o-)ct',^Cj>, (lb) 

with Si = gCTs grCi^s' and sgn(CT) = 1 (—1) for a =t 
(I). Here cj^ creates an electron with spin s =t,4' at 
site i. Si denotes a localized moment of spin S' = 1/2 at 
site i. sgn(i, j) = di x dj/\di x dj\, where di and dj are 
vectors connecting site i and j (see Fig. [T]). Namely, it 
takes 1 (—1) when the electron hops clockwise (counter 
clockwise), respectively. The second term in Eq. ((TJi) 



FIG. 1. (Color Online). Sketch of the lattice. The electron 
of the up-spin state hops from site j to i with itsoSgn(i, j), 
where sgn{i,j) = {di x dj)/\di x dj\. 


denotes the antiferromagnetic (J > 0) coupling between 
the localized spin and the itinerant electron at the same 
site. This model is defined on the honeycomb lattice and 
has two sublattices, A- and B-site. The next nearest 
hopping tao^ whose characteristic form originates from 
the spin-orbit interaction, results in the bulk band gap 
for J = 0, thereby stabilizing a topological insulating 
phase for the conduction electrons. 

In order to analyze the system, we employ DMFT, 
which enables us to treat local correlation exactlyi^^— 
With this approach, we map the system onto an effective 
impurity problem and solve it self-consistently. The self- 
consistency equation is given by 


9a^{^) = X 


1-1 


{U3 + I5)l-K{k)-^R{U3) 


( 2 ) 


where /lo-(k) is the two-dimensional matrix obtained by 
the Fourier transform of the hopping matrix. The self¬ 
energy (^) of the retarded Green’s function for the 
lattice system is a diagonal matrix and can be computed 
from the Green’s function gcr(w) of the effective impurity 


model. In this study, we employ the numerical renormal¬ 
ization group (NRG) method for the computation of the 
self-energy)^^— which is appropriate for calculations in 
the low temperature region. 

For J = 0, our system is reduced to the ordinary Kane- 
Mele model. In this case, the topological structure of the 
ground state is characterized by the spin Ghern number 
.^sCh, which is proportional to the spin-Hall conductiv¬ 
ity. In terms of the Green’s function, this topological 
invariant is defined as follows— 

^sCh = J dwdk sgn{a)tT[G~^{k)ditGa{k) 

G-\k)d,.G„{k)G-\k)dpG^{k)\, ( 3 ) 

where {fj,, v,p = Q, 1, 2) is the anti-symmetric tensor 

taking = x, fe := (fc^,, fc^), 0 := (9,^, Here, 

G„{k) := Gct(*w, k) is the two-dimensional matrix in the 
pseudospin (sublattice) space and defined along the imag¬ 
inary axis. This quantity characterizes the map from 
(w, k) to GL{2, C), and is well-defined even in correlated 
systems.— For a continuous change of the ground state 
wave function, the topological structure changes only 
when the Green’s function does not belong to GL(2,Z). 
Namely, if the topological invariant changes with a con¬ 
tinuous transition, the Green’s function satisfies 

detG“^ = 0, (4a) 

or 

detGa = 0, (4b) 

at a certain point of the space (w, k). Physically Eq. (|4li) 
can be satisfied by gap-closing in the single-particle ex¬ 
citation spectrum. The topological transition satisfying 
Eq. (HJd) is possible only for correlated systems. Without 
correlations, it is equivalent to a diverging band width. 
In the presence of electron correlations, Eq. (H}^) can be 
satished by a divergence of the self-energy. 

Throughout this paper, we set the model parameters to 
t = 1 and tso = 0.1, and take the distance of neighboring 
sites as the unit of length. Eor the Green’s function along 
the imaginary axis, we calculate it for Matsubara fre¬ 
quencies iujn = (2n -I- l)7rr with n S Z and T = O.OOOlt. 
The size of the bulk-gap A is A ~ O.Olt, and thus much 
larger than the temperature, T = 0.0001. 


III. RESULTS 
A. Doniach phase diagram 

Let us start with explaining the Doniach phase dia¬ 
gram. The calculated results are summarized in Eig. [21 
For weak exchange interaction J, we can find an anti¬ 
ferromagnetic topological phase, which is induced by the 
RKKY interaction. This phase changes to the trivial an¬ 
tiferromagnetic phase with increasing J. With increasing 
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FIG. 2. (Color Online). Phase diagram of temperature (T) vs. 
the antiferromagnetic coupling (J). Because of the presence 
of the bulk-gap, a first order transition (solid orange line) 
is observed in the weak-coupling region while it changes to a 
second order (dashed green line) as increasing the coupling. In 
this figure, one can observe the topological phase transition. 
For J < 0.683t a topological antiferromagnetic phase (AFTI) 
emerges. This phase changes into a trivial antifferomagnetic 
phase and an ordinary Kondo insulator as the interaction is 
increased. 


J, the Kondo effect becomes dominant and the topolog¬ 
ically trivial Kondo insulator (i.e., the ordinary Kondo 
insulator) appears. Increasing temperature destroys the 
antiferromagnetic order with a first (second) order tran¬ 
sition for weak (strong) interaction strength J. In the 
following, we discuss the details of the phase diagram. 

First, we focus on the properties of the magnetic phases 
in Fig.[2l If the interaction strength J is weak, the RKKY 
interaction is dominant, and we can find the magnetic 
phase. Besides, if J is sufhciently weak, a first order 
transition is observed as a function of temperature, the 
origin of which can be found in the bulk-gap. The tem¬ 
perature dependence of the magnetization is plotted in 
Fig.|3Ka) and (b) for J = O.lt and J = t, respectively. For 
J = O.lt the magnetization shows a jump at T = O.OlSt 
[see Fig. Ha)]. Here we have defined the antiferromag¬ 
netic moment for the conduction electron and localized 
spin at A-site as me = — ni^i)l2 and ruspin = (5'f) 

with i £ [A-site], respectively. With increasing J, this 
phase transition becomes second order; the Kondo tem¬ 
perature becomes energetically larger than the band gap 
and a continuous transition is observed for J = t [see 
Fig. Hb)]. Correspondingly the quantum phase transi¬ 
tion at Jci ~ 1.61f is also continuous; the magnetization 
around zero temperature {T = O.OOOlf) approaches zero 
continuously [see Fig. He)]. 

Next, we discuss the topological structure of the phase 
diagram. The antiferromagnetic phase can possess topo¬ 
logically nontrivial properties We can use a simpli¬ 
fied formulai^ for the spin Chern number, because the 
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FIG. 3. (Golor Online). The spin Ghern number and mag¬ 
netic moments. In panels (a) and (b), temperature depen¬ 
dence of the antiferromagntic moments is plotted for J = t 
and J = O.lt, respectively, rric and nispin denote the mag¬ 
netic moment of the itinerant electron and localized spin, (c): 
the antiferromagntic moment as functions of the interaction 
strength, (d): the spin Ghern number as a function of the 
interaction strength J. 

Green’s function is not singular in this magnetic phase 
as we will show in the next subsection and the Appendix; 

A^sCh = ^ J dkd^(k) ■ dk^d^ik) x dkyda{k), 

(7 

(5a) 

with 

h'^{k) := ha{k) + = 0), 

= dcr.o(fc)]l + ^ d(7,i(fe)fi, (5b) 

i 

and di = di/\d\. Here 1 and Ti {i = x,y,z) denote the 
two-dimensional identity matrix and the Pauli matrix for 
the sublattice space. We have confirmed that the results 
are consistent with those obtained from the definition 
[Eq. (H]- The results are plotted in Fig. Hd). 

As seen in this figure, we find that the system shows the 
topologically nontrivial properties for J < 0.683t, while 
a further increase of the interaction strength induces a 
topological transition. At this continuous topological 
transition point, gap-closing in the density of states is 
observed; the Green’s function becomes singular, satisfy¬ 
ing Eq. (H-) at iuj = 0. The gap-closing is shown in Eig.S) 
In this figure, the momentum-resolved spectral function 
A{uj,k) ~ —Im^^ trG^(w, A;)/7r is plotted for J = 0.6t, 
0.683t, and 0.8t from bottom to top. The right panels are 
magnified versions of them. As it is observed in the ordi¬ 
nary Kondo lattice model, a renormalized band appears 
around w = 0 because of the correlations between the 
conduction electrons and the localized spins. When ap¬ 
proaching the topological phase transition point from the 
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nontrivial phase (J = 0.6t), we can see that the bulk-gap 
becomes firstly narrower (Jc 2 = 0.683t), and then widens 
again with further increasing of J (J = 0.8t). A dif¬ 
ference in the origin of the gap formation between these 
two phases gives rise to this gap-closing. For J < Jc 2 the 
bulk gap arises from the hopping tso, while for J > Jc 2 , 
the antiferromagnetic order induces the gap. Here, Jc 2 
denotes the critical value for the topological phase transi¬ 
tion. Putting the above results together, we end up with 



momentum k momentum k 


FIG. 4. (Color Online). Momentum-resolved spectral func¬ 
tion A{u),k) = —ImX)o-trC^ (w, fc)/rr for several values of 
the interaction strength J. In this figure the momentum 
k sweeps the high symmetry points in the Brillouin zone; 
{kx,ky) = (0,0) —>■ (7r/%/3, tt) —>• (0,47r/3) —>■ (0,0). Here 
we have taken the distance of neighboring sites as the length 
unit. The data for J = 0.6t, 0.683t, and 0.8t are plotted from 
bottom to top. Right panels are the magnified version around 
cj = 0. 

the Doniach phase diagram shown in Fig. [2] 

We finish this subsection with a comment on the topo¬ 
logical structure in the Kondo lattice. As discussed in the 
next subsection in detail, there is no topological phase if 
the system is paramagnetic because zeros of the single¬ 
particle Green’s function (i.e., pole of the self-energy) 
appear. On the other hand, in Refs. [2^ and [13, the non¬ 
trivial phase persists even in the paramagnetic phase. 
This apparent contradiction arises from the difference of 
the used methods; in Refs. and [13 the localized spin is 
treated with the slave-boson approach and is separated 
into a boson and a fermion which participate in the topo¬ 
logical structure. The authors study the resultant model 
at the mean-field level, obtaining the paramagnetic topo¬ 
logical phase. However, if the local correlation is treated 
exactly, the topological structure in the paramagnetic 
phase is destroyed due to zeros of the Green’s function 
in the limit of the Kondo lattice. On the other hand, 
for small J, the RKKY interaction induces the antifer¬ 
romagnetic order and removes the zeros of the Green’s 
function. Thus, in the Kane-Mele Kondo lattice model, 
the nontrivial state is stabilized only in the antiferromag¬ 


netic phase. 


B. Restoration of topological properties at finite 
temperatures 


In this section we will demonstrate that in heavy- 
fermion systems with topological structure, the inter¬ 
play between electron correlations and topological prop¬ 
erties leads to an intriguing crossover behavior; topolog¬ 
ical properties are restored at finite temperatures in the 
region where the topological structure of the ground state 
is destroyed. This behavior is observed via analysis of 
the bulk and the edge. Calculations for the bulk systems 
reveal that the spin-Hall conductivity increases even if 
the topological structure is absent at zero temperature. 
Calculations under ribbon geometry [the open (periodic) 
boundary condition for x- (y-) direction] reveal that gap¬ 
less edge modes appear due to temperature effects. In 
the following, we will confirm these intriguing properties 
in two steps. 
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FIG. 5. (Color Online), (a): temperature dependence of the 
spin-Hall conductivity, (b) ((c)): self-energy of up-spin state 
in the paramagnetic (antiferromagnetic) phase, respectively. 
Here, we plot them along the imaginary axis. In the antifer¬ 
romagnetic phase, we plot the self-energy at the A-site where 
the electron state for up-spin is majority. 

In order to see the essence of the restoration of topo¬ 
logical properties, let us first restrict ourselves to para¬ 
magnetic solutions. Temperature dependence of the con¬ 
ductivity is plotted in Fig. [SKa). In the low tempera¬ 
ture region, the spin-Hall conductivity is zero since the 
topological invariant is no longer well-defined due to the 
Kondo effect; the singlet formation between electrons and 
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localized spins leads to zeros of the Green’s function (i.e., 
divergence of the self-energy) [see Fig. |5l[b)]42 With in¬ 
creasing temperature, the Kondo effect is suppressed, and 
the conductivity increases. For T > 0.03t, the conduc¬ 
tivity approaches the values of J = 0 which are almost 
quantized. The increase of the spin-FIall conductivity is 
also observed for J = 1.5t, even though increasing the 
coupling strength J suppresses the conductivity at finite 
temperatures due to enhancement of the Kondo effect. 


(a) T = m-H J=\M (b) T/t = 0.1 Jjt = 1.5 



- 0.01 0.00 0.01 




FIG. 6. Momentum-resolved spectral function for up-spin 
state A^(uj,k) = —Im trG:^(cj, fe)/7r at the edge site a; = 0 
for T = O.lt and T = O.OOOlt. In the spectral weight for 
down-spin state, we can find edge modes propagating in op¬ 
posite direction. The R-DMFT is performed under the open 
(periodic) boundary condition for x- (y-) direction. We have 
zig-zag edges at a; = 0 and L—1. For T = O.OOOlt edge modes 
are gapped due to the Kondo effect, while the gap disappears 
with increasing temperature. For T = O.lt, an extremely tiny 
gap is observed. We conclude that this is due to the finite 
size effect because the gap width should be the same as that 
for T = 0.0001 if the Kondo effect induces it. 

This increase of the conductivity is interpreted as a 
restoration of gapless edge modes. Our real-space dy¬ 
namical mean-field theory (R-DMFT) calculations using 
the ribbon geometry reveal how finite temperatures affect 
the edge modes. In Fig. [6] (a) and (c), we can see that 
edge modes are destroyed due to the Kondo effect for 
T = O.OOOlt, which is consistent with the bulk behavior. 
With increasing temperature, the gap in the local density 
of states at the edge site disappears [see Fig. EKc)]. Cor¬ 
respondingly, we can observe a restoration of edge modes 
in Fig. [IKb), which leads to an increase of the spin-Hall 
conductivity. 

We have so far demonstrated the characteristic 
crossover behavior at finite temperatures in heavy- 
fermion systems. The origin of this crossover is the com¬ 
petition of the Kondo effect and the topological prop¬ 
erties. In the low temperature region, the Kondo effect 
governs the low energy properties and destroys the topo- 


J/f=-0.5 J/f=-1.5 ■ 

J/f=-0.75 ^ d/f=-2.0 
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FIG. 7. (Golor Online). Temperature dependence of the spin- 
Hall conductivity for several values of coupling strength J. 
The error bars in the low temperature region arise from ex¬ 
trapolation to the oj —>■ 0 limit. At high temperatures these 
error bars are smaller than symbols. Arrows denote antifer¬ 
romagnetic transition points. 

logical structure. On the other hand, topological proper¬ 
ties (i.e., the conductivity and edge states) are restored 
if the temperature is higher than the Kondo temperature 
but smaller than the energy scale of the band gap of the 
non-interacting topological insulator. 

An important point to be noticed here is that the 
restoration of topological properties is also observed in 
the antiferromagnetic phase. Namely, the increase of the 
spin-Hall conductivity is also observed in the antiferro¬ 
magnetic phase. In Fig. [71 temperature dependence of 
the conductivity in the antiferromagnetic phase is plot¬ 
ted. For J = 0.75t, the ground state is topologically 
trivial. Correspondingly, the spin-Hall conductivity van¬ 
ishes at zero temperature, but increases with increasing 
temperature. We note that the ground state is an antifer¬ 
romagnetic topological phase for J < 0.683. As seen in 
Fig.EKc), the magnetic order removes the pole of the self¬ 
energy, and the ground state possesses nontrivial prop¬ 
erties (properties of the self-energy in the trivial anti¬ 
ferromagnetic phase are discussed in Appendix [^. This 
topological phase leads to a dip structure in the tem¬ 
perature dependence of the conductivity for J = 0.5t; 
with decreasing temperature the Kondo effect firstly de¬ 
creases the conductivity for 0.02It < T < 0.03t, but with 
entering the antiferromagnetic phase (T < 0.02It), the 
conductivity increases again because the ground state is 
the antiferromagnetic topological insulator. 


IV. CONCLUSION 

In this paper, we have analyzed the Kane-Mele Kondo 
lattice model at finite temperatures and obtained the Do- 
niach phase diagram under the topologically nontrivial 
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condition (Fig. [2]). 

In particular, we have observed an intriguing crossover 
behavior due to the interplay between electron correla¬ 
tions and topological properties. In the low temperature 
region, the Kondo effect destroys the topological struc¬ 
ture, while the topological properties are restored if the 
temperature is higher than the Kondo temperature but 
smaller than the band gap of the topological insulator 
in the non-interacting case. The restoration of topo¬ 
logical properties can be observed in the bulk and the 
edge. The spin-Hall conductivity rapidly increases even 
if it is almost zero in the low temperature region. Edge 
modes, destroyed by the Kondo effect in the low temper¬ 
ature region, appear with increasing temperatures. The 
crossover behavior is observed in both of the paramag¬ 
netic and the antiferromagnetic phases. 
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Appendix A: SELF-ENERGY IN THE 
ANTIFERROMAGNETIC PHASE 

In this appendix, we discuss why the topological anti¬ 
ferromagnetic phase is not stabilized around the quan¬ 
tum critical point (Jci ^ I.61t) in spite of the small 
magnetization. By using the Hartree approximation, the 
self-energy is estimated as E-j- ^ Jiric in the antiferro¬ 
magnetic phase. Furthermore, around the quantum crit¬ 
ical point Jci ^ 1.6It, the magnetic moment approaches 
zero. With these facts, one might expect that a topolog¬ 
ical phase can be stabilized. This is, however, not the 
case. Although the self-energy is non-singular in the an¬ 
tiferromagnetic phase, a huge real-part of the self-energy, 
—ReE-|-, drives the systems into the trivial phase around 
the critical point. In the following we will discuss details. 


First, in the absence of a singularity of the Green’s 
function the topological structure is well-defined. In 
Fig. (HKa), the self-energy for J = 0.8t is plotted. In 
Figs. ISKc) and UKa) the self-energy shows no divergence 
in both the topologically nontrivial and trivial phases. 
With increasing J, —ReE-|'(ia;o) becomes larger, which 
induces the gap-closing in the antiferromagnetic phase 
at Jc 2 ~ 0.683t. 

Second we show that the real part of the self-energy, 
—ReE-|-(*u;o)j rapidly increases when approaching the 
quantum critical point Jd, see Fig. |S] (b). This figure 
indicates that — ReE-f(*a;o) increases rapidly with ap¬ 
proaching the critical point even though the magnetic 
moment becomes small. The large value of —ReS-i'(iaJo) 
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FIG. 8. (a) self-energy at J = 0.8t. (b) the real part of the 
self-energy at iujo. 


induces a bulk-gap and destroys the gap induced by the 
next-nearest-neighbor hopping at Jc 2 = 0.683t. Hence, 
we cannot hnd a topological phase around the quantum 
critical point. 

Therefore, we conclude that the system is in a topo¬ 
logically trivial phase around Jd because of a huge value 
of —ReE(za;o). 
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